Critical Behavior of CP^ at 6* = vr, Haldane's Conjecture and the Universality Class 
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Using an approach to analyze the 8 dependence of systems with a ^-term we recently proposed 
, the critical behavior of CP^ at S = tt is studied. We find a region outside the strong coupling 
regime where Haldane's conjecture is verified. The critical line however does not belong to the 
universality class of the Wess-Zumino-Novikov-Witten model at topological coupling k — 1 since it 
shows continuously varying critical exponents. 
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Two dimensional CP^^^ models are of great interest 
in high energy physics. They are toy models for testing 
nonperturbative and topological properties in a confin- 
ing, asymptotically free quantum field theory and share 
also with QCD, the theory of the strong interaction be- 
tween particles, a well defined topological charge, possess 
instanton solutions, have 0- vacua and admit a 1/N ex- 
pansion. Understanding the role of the 6 parameter in 
QCD and its connection with the "strong CP problem" is 
one of the major challenges for high energy theorists 1]. 
Unfortunately, Euclidean lattice gauge theory, our main 
nonperturbative tool for QCD studies, has not been able 
to help us for two reasons: (i) the difficulties in defining 
the topological charge operator on the lattice Q and (ii) 
the imaginary contribution to the action coming from the 
6 term that prevents the applicability of the importance 
sampling method. CP^~^ models regularized on a lat- 
tice have a well defined topological charge but share with 
QCD point (ii) and this means that, from the point of 
view of numerical calculations, the complex action prob- 
lem is as severe as in QCD. Indeed most properties of 
QpN-i T^Q^^Qig have been obtained in the context of the 
1 /N expansion 0, |3| , showing a qualitative behavior as 
a function of similar to that of the massive Schwinger 
model at weak coupling Q. The cusp in the free energy 
density at 6 = tt signals a first order phase transition 
at this value of 9 with spontaneous breaking of CP, in 
agreement also with the strong coupling results Q . This 
behavior, which is expected to persist for finite values 
of N, at least until = 4, is common to most of the 
models the 9 dependence of which is known, and was 
conjectured for non-abelian gauge theories by 't Hooft 
in 1981. Some time ago, using different arguments 
and under very general assumptions, we also argued that 
a singular behavior in 9 is to be expected in systems with 
a 0-term 0]. For A^ < 3 dislocations break the asymp- 
totic scaling (continuum limit) Q and the 9 dependence 
becomes more intriguing in these cases. 

However high energy physics is not the only field where 
topological structures and the 9 term play a fundamental 



role. In condensed matter physics Haldane showed [lO| 
that chains of quantum spins with antiferromagnetic in- 
teractions, in the semiclassical limit of large but finite S, 
are related to the two-dimensional 0(3) nonlinear sigma 
model at coupling g = 2/[S{S + 1)]^^^ and topological 
term 9 — O^n (integer, half-integer spin). While integer 
spin chains have a mass gap, half-integer chains should 
be gapless. Haldane conjectured that the 0(3) non linear 
sigma model presents a second order phase transition at 
9 = 7:, keeping its ground state CP symmetric. Based 
on a partial summation of the strong coupling expan- 
sion, Affieck ll| argued that the two-dimensional 0(3) 
non linear sigma model has indeed a gapless phase at 
9 = TT and weak coupling. Furthermore Affleck and Hal- 
dane fl5| also argued that the critical theory for generic 
half-integer spin antiferromagnets is the Wess-Zumino- 
Novikov-Witten model at topological coupling fc = 1, 
with the following scaling law for the mass gap |13i] near 
the phase transition point: 



m{9) 



{Tr~9f^^ln[{Tr-9)-^/'^] 



(1) 



In 1995 Bietenholz, Pochinsky and Wiese [IJ] per- 
formed a simulation of the two-dimensional 0(3) non 
linear sigma model with a 0-term in order to verify Hal- 
dane's conjecture. They used an efficient Wolff cluster 
algorithm [l5| on a triangular lattice with a constraint 
in the action which preserves the 0(3) symmetry, and 
simulated the system at a fixed value of the coupling 
(temperature) g = oo which, due to the constraint in 
the action, should be outside the strong coupling regime. 
By measuring the topological and magnetic susceptibil- 
ities and using finite size scaling theory, the authors of 
[li^ found a second order phase transition at 9 = tt con- 
firming Haldane's conjecture, and a finite size scaling in 
relatively good agreement with the assumption that the 
critical exponents of the phase transition are those of the 
k = 1 Wess-Zumino-Novikov-Witten model ([T]). 

As stated before, the sign problem has restricted very 
much the lattice non-perturbative investigations of sys- 
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terns with a 9-teim. A few years ago we introduced two 
alternative schemes 3 121 to analyze the 9 dependence 
in these systems. The first approach [3l was based on 
a precise determination of the probability distribution 
function of the topological charge whereas the second one 
17| . originally inspired by analytical properties of the 
Ising model with an imaginary magnetic field, is very ef- 
ficient to determine the critical behavior at = tt in cases 
in which a continuous transition at this value of 9 shows 
up. The two approaches were successfully tested in sev- 
eral integrable models and, after this consistency checks, 
applied to the analysis of the 9 dependence of the CP^ 
model. Indeed reference [3] contains the first full re- 
construction of the continuum 9 dependence of a CP^^^ 
model, showing in particular the spontaneous CP sym- 
metry breaking at = tt in the continuum limit. Here 
we want to use the second of the previously cited ap- 
proaches 17[ in order to analyze the phase structure in 
9 of CP"'^ which, as well known, is equivalent to the non 
linear 0(3) sigma model. For this reason in the following 
we will briefly summarize the main algorithm steps. 

We have adopted for the action the standard "auxiliary 
U(l) field" formulation 

Sg — 2/^ ^ ^\ ^n-\-fi^nUn.ii ~1~ ^n^n+^Uyi^^ 2) (2) 

where z„ is a 2-component complex scalar field that sat- 
isfies z„z„ = 1 and Un,^^ is a U{1) "gauge field". The 
topological charge operator is defined directly from the 
C/(l) field: 



2ti 



(3) 



where Up is the product of the U{1) field around the 
plaquette and — vr < \og[Up) < tt. 

Our numerical approach uses as input the results from 
numerical simulations of CP^ at imaginary values oi 9 ~ 
~ih (h real) [iTj. The action to simulate S = Sg + 
Se=-ih is then real and local. This implies that the main 
computational cost is practically equivalent to standard 
simulations of CP^ model at 9 — 0. This is important 
since it allows to perform large volume simulations. 

We use x{h), the topological charge density, and the 
variable z — cosh-| to define the quantity y{z) = 
x(h)/ tanh^. Using the transformation y\{z) = y{ei z) 
and plotting y\/y against y one gets typically a smooth 
function for small y [13] ■ Hence we can expect a simple 
extrapolation to y — > {i.e. in the region corresponding 
to real 9) to be reliable and thus obtain x{9) (of course the 
result has to be independent of the specific value of A) . In 
the same spirit the effective exponent 7a = j\og{y\/y) 
as a function of y will give the dominant power of y{z) as 
a function of z near z = or, equivalently, the behavior 
of x{9) for 9 TT. If 7A(y ^ 0) = 1 CP symmetry is 
spontaneously broken at 9 — tt, values between 1 and 2 



FIG. 1: Order parameter in Ising model (full symbols) and a 
toy model discussed in the text (open symbols); the lines are 
analytical results 



indicate a second order phase transition with a divergent 
susceptibility and so on 



17[. Our approach assumes im- 



plicitly that the theory has not phase transitions at 9 < tt 
and in particular that CP symmetry is realized at = 0, 
the last being a necessary condition for the theory to be 
well defined at 6* 7^ 0. 

Fig. 1 shows the goodness of our method when ap- 
plied to two models, the analytical solution of which is 
known [l7| . The first one is the one-dimensional Ising 
model within an imaginary magnetic field, which breaks 
spontaneously CP at 9 — tt. The second one is a toy 
model with only one effective degree of freedom. This 
model, which resembles very much the free instanton gas 
model, has the partition function Zv{9) = (1 -I- Acos9)^ 
for V degrees of freedom. We have plotted in Fig. 1 the 
order parameter x{9) against 9 for both models. The 
lines stand for the exact analytical solution whereas the 
symbols correspond to the values extracted from our ap- 
proach. 

Let us resume our results for CP^ . This model is equiv- 
alent, at /3 = (strong coupling), to two-dimensional 
compact QED with topological charge. CP symmetry 
is therefore spontaneously broken at 9 — tt m this limit. 
This result remains qualitatively unchanged when mov- 
ing to larger (3 values, until /3 ^ 0.5. For (3 > 0.5 CP 
symmetry is recovered at 9 = tt and this result seems to 
hold for any /3. For (3 values between 0.5 and 1.5 the 
order parameter x{9) vanishes at 9 = tt as {tt — ffy^^") 
with e{f3) varying continuously between and 1. This 
implies a divergent topological susceptibility xt = 
at 9 — TT for 0.5 < /? < 1.5 and we get therefore a line of 
second order phase transition points with continuously 
varying critical exponents. For [3 > 1.5 we find e(/3) = 1. 
Numerical evidence of these results is given in figures 2-5. 

The details of the simulations are: we used lattices 40^ 
for (3 < 1.0, 100^ for 1.0 < /3 < 1.5 and 200^ for /? > 1.5 
checking absence of finite volume effects; this choice has 
been taken looking at existing data for the correlation 
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FIG. 2: Order parameter for CP^ model; /3 = 1.0 , L=100. 
Statistical errors are of the order of the symbol size. 



FIG. 4: Effective exponent 7a for /3 = 1.0, L=100 (circles) 
and P = 1.6 , L=200 (crosses) and two different extrapolations 
of /? = 1.0 data. 
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FIG. 3: The same data of Fig. 2 in log-log plot; the line is 
the fit with the funcion const + a[-n — Oy'^^'^ 

lenght at 6* = [IB] (C < 3 for /3 < 1.0, 3 < ^ < 20 for 
1.0 < j3 < 1.5). For each value of /3 we have performed 
simulations for 25-40 values of h with statistics exceeding 
500000 MC iterations each in order to evaluate x{h). 

In Fig. 2 we report the order parameter as a function of 
6 for /? — 1.0, hence in the region where < e < 1; in Fig. 
3 the same data are plotted in a log-log scale and, to give 
the feeling of the accuracy in the determination of the 
critical exponent, a fit with the function x{6) — const + 
a(7r — Oy^f^") is also shown. The fitted parameters are: 
const=0, a=0.441(l) and e(l.O) = 0.736(1). Equivalent 
graphs have been obtained for the other values of f3 in 
the range 0.6 - 1.5 while, for f3 < 0.5, we have a different 
behaviour of the order parameter, approaching a non zero 
value at 9 — IT. For the two larger values of /3, namely 
1.5 and 1.6, we obtain e — 1 and an order parameter 
indistinguishable from a sin(0) function between 9 = 
and 9 ^ n. 

Starting from the same raw data (the simulations at 
imaginary 9) we can do a consistency check performing 
a different analysis using the effective exponent 7a. We 
are interested in the effective exponent in the y — > 
limit and, in Fig. 4, we can see the 7a extrapolation for 



(3 = 1.0 {L = 100 data, lower curve) and for (3 = 1.6 
(L = 200 data, upper curve). For the [3 = 1.0 case 
extrapolations obtained using two different fitting func- 
tions, namely j\{y) = exp(a + by + cy^) and a + by + cy^ 
are shown; the first function is able to slightly better 
reproduce the data in the full P range with respect to 
the simple polynomyal form. The small difference in the 
extrapolated values does not change qualitatively the re- 
sult. 

For the (3—1.6 data there is essentially no need for 
extrapolation, being very clear that the results point to 
the value 7 = 2 in the y = limit. Data shown in 
the figure are those relative to A = 0.5 but exactly the 
same conclusions can be reached using other values of A. 
It is important to notice how, moving to the region of 
smaller (3, the same fitting functions continue to describe 
very well the data and have been used to extrapolate the 
effective exponent at zero y. 

Fig. 5 resumes all results for exponents e and 7 as a 
function of /?: it is clear from this figure that, considering 
the fate of CP symmetry at 6* = tt, the system passes 
from a phase in which the symmetry is broken (/? < 0.5) 
to a symmetric phase (0.5 < (3 < 1.5) where a second 
order phase transition with varing critical exponents is 
present, ending with a region in which CP symmetry is 
still satisfied but without a divergent correlation lenght. 

The location of the point where the system becomes 
CP symmetric &i 9 = -k can not be obtained from our 
data with a high accuracy since different analysis gives 
slightly different results but the existence of such a point 
is evident from both analysis presented here. 

Nevertheless, in our opinion, the important physical 
fact is the existence of a line of second order transition 
points with continuously varing critical exponents. Using 
the hyperscaling hypothesis we can relate e(/3) with the 
critical exponent v for the mass gap. Indeed, following 
Kadanoff, one can assume that the singular part of the 
free energy density near the phase transition points is 



4 



□ 



0.2 0.4 0.6 0.8 1 1.2 1.' 



FIG. 5: Critical exponent e (squares) from fits to the order 
parameter behaviour and effective exponent 7a (actually 7a — 
I in order to be directly comparable with e , triangles). 



proportional to S, being the inverse mass gap, and 
this implies that m{9) ^ [tt — Oy . As a consequence 
simple algebra gives the relation e = 2v — 1. 

In conclusion we have analyzed the critical behavior 
of CP^ at 6 = TT and found a region in the coupling P 
where Haldane's conjecture is verified. The critical line 
shows continuously varying critical exponents and does 
not belong therefore to the universality class of the Wess- 
Zumino-Novikov-Witten model at topological coupling 
k = 1. From the point of view of Quantum Field Theory 
these new critical points open the possibility for new non- 
gaussian fixed points with anomalous dimensions, where 
a non trivial continuum limit could be defined. 
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